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I. 


INTRODUCTION 


The  problem  of  scattering  by  wedges  has  been  treated  extensively 
by  many  authors  [1,2].  Usually  the  analysis  involves  asymptotic  methods 
to  evaluate  far-fields,  with  little  emphasis  on  near-field  calculations. 

It  is  the  near-field  aspect  with  which  the  present  paper  is  concerned. 

An  interest  in  the  near  zone  fields  of  a source  near  an  edge  arises 
when  formulating  a moment  method  solution  to  the  problem  of  a wire  antenna 
mounted  near  the  edge  of  a plate  or  wedge.  In  a previous  paper  [3],  the 
authors  presented  a moment  method  technique  whereby  problems  with  wires 
mounted  .IX  or  greater  from  an  edge  may  be  treated.  Desiring  to  extend 
this  technique  by  synthesizing  an  expansion  mode  to  treat  the  case  of 
the  wire  near  an  edge,  the  analysis  presented  here  forms  a first  step. 

In  the  present  analysis  an  electric  line  source  parallel  to  an 
edge,  with  arbitrary  current  flow,  is  considered.  Simple  expressions  for 
the  near  zone  fields  are  derived  which  exhibit  not  only  the  correct  edge 
behavior,  but  the  source  singularity  as  well. 


II.  THEORY  FOR  PERFECTLY  CONDUCTING  WEDGE 

A.  Introduction 

The  geometry  of  the  problem  is  shown  in  Figure  1.  The  field  point 
is  (p , 4>)  while  the  electric  line  source  is  at  (p',<t>').  The  line  source 
is  parallel  to  the  edge  of  the  wedge  and  has  components  Iq,  1$,  and  1^ 
of  current  in  the  p,  <j>,  and  z directions,  respectively.  The  current  is 
assumed  to  be  constant  with  respect  to  z.  Thus  the  source  current  density 
may  be  written  as 

T(p,*)  - (IpP+y+Izz)  i 6(p-p')6U-#')  A/M2.  (1) 

The  method  of  solution  is  to  use  the  two-dimensional  dyadic  Green's 
function  for  the  wedge  to  find  the  magnetic  field  due  to  the  vector 
current  T.  Then,  since  the  field  and  source  points  are  close  to  the  edge, 
small  argument  approximations  are  used  for  the  Bessel  functions  which 
occur  in  the  resulting  infinite  series.  The  series  is  then  summed  in 
closed  form,  preserving  the  source  singularity.  Maxwell's  equations  are 
used  to  find  the  F fields.  Of  course,  with  the  fields  of  an  electric 
source  known,  the  fields  of  a magnetic  source  may  be  found  from  duality. 

B.  Two-Dimensional  Dyadic  Green's  Function 

The  two-dimensional  dyadic  Green's  function  is  a solution  of  the 
vector  wave  equation 

(vxvx-k2)6  9(R IR-1  )=vx[!6(R-R')]  (2) 
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Figure  1.  Geometry  of  the  line  source  near  an  edge. 

o _ ^ ^ _ 

where  k = y~  , I is  the  unit  dyad,  !T=pp+<t>4>,  TT*  =p  'p+<p'^,  and  5^  satisfies 
o ^ 

the  boundary  condition  nxVxG^O  on  s.  The  fields  are  then  found  from 


H(R)  = J ^(RlW)  * T(1T')  ds,  (3) 

and,  assuming  eJayt  time  dependence, 

ron  = iivxmn)  - vmi  . (4) 

j cot  ^ 


The  integration  in  Equation  (3)  is  over  the  source  region.  Using 
Equation  (1)  in  Equation  (3)  gives 

m » ip  5m2(*|in-p  + i^m2(ir|ir').^izGm2(Tr|ir').z  . (5) 

The  Green's  function  Gm2  is  given  by  Tai  [4,5]*: 


*Note  that,  from  Equation  (1),  v*Gm2=0  so  that  only  the  FT  and_N  vector 
wave  functions  are  needed  in  the  eigenfunction  expansion  of  Gm2. 
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6m?(R|in  = 


U l 

) n= 


where  6q  = jj  J'j}  , c=nv=n  j^-y 


- ff  , TT ' + R . ff' 

r OaA  oaA  eaA  eaA 


(1+60)(x  -k  ) 


; n=0,l,2***  and  the  vector 


wave  functions  are* 


(Ap)  "la  ♦£ 


Ne  , = ♦ z 


A J'CApJrJJa 
ot  'sin 


Consistent  with  Tai's  notation  [4],  FT'  and  N'  in  Equation  (6) 
refer  to  FT  and  IT  with  primed  spatial  (source)  coordinates,  whereas 
elsewhere  the  prime  refers  to  differentiation  with  respect  to  the 
argument. 

In  defining  field  components,  superscripts  refer  to  the  source 
polarization  and  subscripts  refer  to  the  field  component.  Then,  from 

Equation  (5-7),  H*  = Hp  = H*  = = H*  = 0,  and 

p P <t>  <t>  z 


H^(p,4>)  = TFT 


tt-<(>0;p 


00 

u 1 

* n= 


">  aAJ  (Ap)j  (Ap1) 

T 2 2 — _ — cosoi^si  na4> ' (8) 

=0  (l+6  )(A  -k  ) 


* -2t  r - a20  (Ap)J'(Ap') 

(p  * 4* ) = IV  ^ ."V  1 dA  £ k — w—  cosa<f>cosa<j>'  (9) 

Z (2*“V  ) n=0  (l+60)(AZ-kZ) 


*The  2-D  form  of  the  vector  wave  functions  given  above  is  most  easily 
found  from  the  3-D  form  given  in  Reference  [4]  by  letting  h=0. 
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Hp(p,<f>)  = 


QQ 

[ r 

dx  l 

^\b  j ni 


C “ aXj  (Xp ) J (Xp  1 ) 

\ dx  j 2 2_ — COSa<)>Si na<}i 1 

3 n=0  (l+60)(xZ-k2) 


-21  r « X2J'(xP)J  (XP‘) 

H:(p.<I>)  = vrK  t \ dX  7 2_ — sina^sina^ ’ . (11) 

♦ ) J n£0  (1+5  )(X2-k2) 

* 0 


The  integrations  in  Equations  (8-10)  can  be  evaluated  [4]  to  give 


, I co  Ja(kp)H^(kp')  p<p* 

,,P/  _ J P V aCOSa<j)Sina<(.'  ) 

HZ(p’^  " (2tt-(J)  Tp'  L n+6  T"  \ /o\ 

0 n'°  0 H^2)(kp)Ja(kp’)  p>p' 


H*(p,*)  = 


V COSa(j)COSa()) ' 


V n=0 


„Z,  Z v aCOSa()>Si  na<j>' 

Mp»*J  - n-A  -in  A. n+x  i 


:ir'Vp  n=0 


I|Z/  _ JirkIz  “ sina<j>sina<f>' 

V"'*’  - TS^TT  n£0  TT^r- 


J0(k»)«<21'(kp’)  p<»' 


Ha  (kp)°a(kp,)  P>P’ 


Ja(kp)Hi2)(kp,)  P<P' 


H^;(kp)Ja(kp')  p>p' 


j;(kP)H^;(kp')  p<p' 
(kp)J  ( kp  ' ) p>p‘ 

a a 


Up  to  this  point  the  field  expressions  are  completely  rigorous. 
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C.  Near  Field  Calculation 


The  source  and  field  points  are  now  restricted  to  be  close  to  the 
edge  of  the  wedge,  so  kp <1  and  kp ' < 1 . Using  small  argument  approximations 
for  the  Bessel  functions  in  Equations  (12-15)  and  trigonometric  identities 
yields 


Hp(p’8)  ' (2V-b» 


H?(p,»)  • z 


(27T-4>o)p 


/ °° 

7 I (rr)a[si na( 4>+<j>' ) -sin a( <(>-<}> 1 ) ] 
c n=l  p 

< 

oo 

7 1 — J^Csi na( d>+<#> ' )-sina(4>-4>* 

n=l  p 


7 I (~r)a[c0Sa( 0+$'  )+cosa(<f>-<)> 1 )] 
c n=0  p 

< 

oo 

~ i-  I (~ — ) aC cosa ( 4)+<J>  1 ) +cos a ( <f>— ' )]+l 
V.  c n=0  p 


i 08 

7 l (~r)a[sina(  ((>+(() 1 ) -si na ( <j>-4> ' )] 
c n=l  p 

I 7 I (~ ) si  not  ( 4>+<i> 1 )-sina(  <(>-  <(> 1 )] 

n=l  p 


( i * 

7 l (^-r)°t[cosa(  4>+<t> ' ) -cos a ( 4>- 4> ' ) ] 
n=l  p 


[-7  l ) a[cosoi( <|>+4> 1 )-C0Sa(<J>-4>'  )] 
w n=l  p 


p<p ' 

p>p ' 

(16) 

p<p  ' 

P>P  1 

(17) 

p<p ' 

p>p' 

(18) 

p<p ' 

p>p' 

(19) 
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Using  the  summation  formulas  A3  and  A4  of  the  Appendix  gives,  after 
simplification,  the  desired  closed  form  results: 


The  If  fields  could  be  found  either  by  using  (4)  on  (12-15)  and 
applying  the  procedure  of  small  argument  approximations  and  summation 
of  the  series,  or  using  (4)  on  (20-23)  directly.  Identical  results 
are  obtained  with  either  method  for  all  F components  except  for  E|, 
which  would  reduce  to  zero  if  the  curl  of  (22,23)  were  taken.  Thus, 

Ey  is  found  from  (4)  and  (14-15),  while  the  other  E components  are 

found  using  (4)  and  (20,21).  After  simplification,  E^  = E£  = E‘  = E|  = 0, 

and 


W A 


I v(pp')V*1 

E*(p,*)  = 2j^^v! 


(p  '^v-p^v)sinv(<J>+4>' 

L < 


(P|2v-p2v)sinvU-*') 

n2 


E*(P,*) 


i^v(pp')V~1 

2 j toe  Q ( 2 7T  - <J>Q ) 


/ 2v  ,2vv 


p p 

siny(<}>+4> 1 ) + (p  v-p  1 ^sinvU-ji1 ) 


E^(p  ,$) 


. V(pp')V~  (2ovp'v-(p'2%2v)cos,[«tf) 


2 j oie  0 V 2 tt—  4>0  J 


+ 2pvp  |V-(p  |2v+p2v)cosvU-<)>1  )~]  _ _j 


J(%p 


6 (p-P  ' ) 6 ( 4>- <J> 1 ) 


From  Equations  (4,14,15), 


z -*k2lz 

Ez(p^}  = ueo(yiT-voy  ^ sino*sin«*’ 


J (kp)H*2)(kp') 


a a 


JweQp 


i(p-p  ' )6  (<>>-<}) 1 ) 


ot  (kp )^a(kp  ' ) p>p' 


Using  small  argument  approximations  in  the  Bessel  functions  in 
Equation  (28)  gives 


z A I"1 

Ez‘p"»  ■ 


l C cos ct ( 4>+<#> 1 ) —co Sc* ( 4>—  <f> ' )]  ~ (2t-)“ 


Jw£op 


> (p -P  ' )6  ( 4>-<f> 1 ) 


l C COS  a ( <J>+4> 1 ) -COSa(  <)>-<(> 1 ) ] 1 (M 

n=l  n p 
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Using  the  summation  formula  (A7)  of  the  Appendix  gives 


6(p-p  ’ )6(<p-<t>') 


(30) 


Equations  (20-23),  (24-27),  and  (30)  are  the  complete  field  expressions 
for  an  arbitrary  vector  electric  line  source,  and  are  rigorous  in  the 
near  field  limit,  including  the  proper  source  singularity.  Note  the 
proper  edge  behavior  in  the  fields  as  p+O.  Also  note  that  the  expressions 
for  H in  Equations  (20-23)  are  independent  of  frequency,  indicating  a 
quasi-static  result.  Hence  (20-23)  could  have  been  obtained  from  a 
conformal  mapping  approach. 

Figure  2 shows  the  magnitude  of  E^  vs.  distance  along  the  surface 
of  the  wedge  for  various  wedge  angles  and  x=lm,  p'  = .lm,  4> ' =0 , and  1^=1. 

Although  the  transverse  field  components  are  singular  at  the  edge 
of  the  ideal  wedge,  no  such  singularity  exists  on  any  physically 
realizable  electromagnetic  scatterer.  The  reason  for  this  difference  is 
that  the  ideal  wedge  has  a perfectly  sharp  edge  with  a truly  discontinuous 
normal  vector,  while  in  physical  reality  any  real  "edge"  will  have  a 
small  radius  of  curvature  so  that  the  normal  vector  at  the  "edge"  will 
change  continuously  as  it  moves  around  the  "edge".  This  fact  was 
recently  used  by  Rhoads  [6]  in  connection  with  physically  realizable 
antenna  aperture  distributions. 

Another  interesting  observation  is  that  some  field  components 
become  singular  when  the  line  source  approaches  the  edge  (p'-K)).  This  is 
in  agreement  with  reciprocity.  Of  course,  this  situation  could  never 
happen  in  practice  because  the  perfectly  sharp  wedge  does  not  exist. 

In  order  to  get  a feeling  for  the  field  behavior  at  a physical 
"edge"  the  mathematically  ideal  edge  of  the  wedge  of  this  section  is 
removed  by  adding  a circular  cylinder  tip.  Although  this  geometry  is 
not  an  exact  model  of  a physical  edge,  it  should  allow  a qualitative 
view  of  the  field  behavior,  with  the  advantage  that  it  is  mathematically 
tractable. 


III.  THE  CYLINDER  TIPPED  HALF-PLANE 

In  this  section,  using  the  same  method  as  outlined  in  Section  II, 
the  closed-form  field  equations  for  a cylinder  tipped  wedge  are  treated. 
Only  the  expression  for  E|  for  the  special  case  of  a half-plane  (<f>o=0) 
with  a I ^ line  source  (Iz=Ip=0)  on  the  surface  of  the  half-plane  (<j>'=0) 
is  presented  here,  but  the  general  case,  for  all  field  components  and 
arbitrary  wedge  angle,  could  be  developed  in  the  same  manner. 
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The  geometry  for  the  cylinder  tipped  half-plane  is  shown  in 
Figure  3.  The  two-dimensional  dyadic  Green's  function  may  be  obtained 


Figure  3.  Geometry  of  an  electric  line  source  on  the  surface  and 
close  to  the  edge  of  a cylinder^tippeg  half-plane. 
Current  polarization  is  in  the  4 (or  y)  direction. 


from  the  3-D  result  given  by  Tai  [4]  as  described  in  Section  II.  The 
Green's  function  is  seen  to  consist  of  a sum  of  the  half-plane  Green's 
function  and  terms  which  account  for  scattering  from  the  cylindrical 
tip.  Thus, 

r»(R)  = r»(wed<ie)(R)  + 


where  p*»(wed9e)(]f)  is  given  by  Equations  (26,27).  Omitting  the  detailed 
derivation. 


The  total  Ej  field  for  y=0,  *'=0,  $Q=0,  p'  = .lm,  A=lm  is  plotted  vs. 
p=x  for  various  values  of  a for  the  cylinder-tipped  half-plane  in  Figure 
4.  Also  shown  in  this  figure  is  the  ET  field  for  a perfectly  sharp  half- 
plane, which  would  correspond  to  a=0  (rto  cylinder  tip).  Figure  i shows  the 
same  curves  for  p '= .01  (A=lm).  The  p=x  axis  is  logarithmic  to  give  an 
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Figure  5.  Magnitude  of  e!  (phase  = -90°)  versus  position  (p=x)  along 
the  surface  of  a cylinder  tipped  half-plane  for  various 
values  of  cylinder  tip  radius,  a.  p'^.Olm,  A=lm,  <(>'=0, 
<j>=0,  I =1,  <p  =0  (half-plane). 


expanded  view  of  the  edge  area.  First  of  all,  note  that  for  a=0  (no 
cylinder)  the  field  is  becoming  singular  as  expected.  Now,  for  a small 
but  finite  cylinder  radius  a,  the  field  increases  as  the  edge  is 
approached,  but  reaches  a maximum  and  goes  to  zero  at  p=a,  as  required 
by  boundary  conditions.  It  can  be  seen  from  Figures  4 and  5 that  the 
field  reaches  a higher  maximum  for  smaller  a,  but  is  not  singular  except 
for  a=0. 

Figure  6 shows  the  total  e|  field  for  y=0,  <j>'=0,  4>o=0,  p'  = .lm, 

A=lm  and  a=5xl0-^m  vs.  p=x  on  aVectangular  scale  for  tne  ideal  half-plane 
and  the  cyl inder-tipped  half-plane,  along  with  a sketch  of  the  relative 
size  of  the  cylinder.  Figure  7 shows  the  same  curves  for  a=5xl0"5m. 


IV.  CONCLUSION 

Simple  closed  form  expressions  have  been  derived  and  presented  for 
an  electric  line  source  with  arbitrary  current  flow  near  and  parallel 
to  the  edge  of  a perfectly  conducting  wedge. 

The  effect  of  a physical  edge  (not  perfectly  sharp)  has  been 
investigated  by  comparing  the  fields  near  the  end  of  a cylinder-tipped 
half-plane  with  those  near  the  edge  of  an  ideal  half-plane. 

« 

The  above  results  indicate  that  incorporating  the  singular  edge 
behavior  into  a solution  may  not  be  required  to  obtain  agreement  with 
experiment.  However,  if  one  is  solving  an  integral  equation  based  upon 
the  ideal  wedge,  incorporating  the  singular  edge  behavior  may  improve 
the  stability  and  convergence  of  the  solution  [7], 


CYLINDER 

'TIPPED 


Figure  6.  Magnitude  of  (phase  = -90°)  versus  position  (p=x)  along 
the  surface  of  a cylinder  tipped  half-plane  compared  with 
that  for  an  ideal  half-plane.  The  relative  size  of  the 
cylindrical  tip  is  also  shown.  p'=.lm,  X=lm,  <f>'=0,  d>=0, 

I =1,  a=p '/200. 
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Figure  7. 


Magnitude  of  (phase  = -90°)  versus  position  (p=x)  along 
the  surface  of  a cylinder  tipped  half-plane  compared  with 
that  for  an  ideal  half-plane.  The  relative  size  of  the 
cylindrical  tip  is  also  shown.  p'=.lm,  A=lm,  <f>'=0,  $=0, 

I =1,  a=p  72000. 
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APPENDIX 

SUMMATION  OF  SERIES 


Starting  with  the  known  result 


V anx  _ J 

n=0  1-a 


replace  a by  aejl*’: 


x>0,  0<a<l  , 


V nx  jn<i>x  _ 1 _ l-axe”^Xl,> 

n=0  ‘ l-aV*x  ' l-2axcos(x*)«Zx  ‘ 

Equating  real  and  imaginary  parts  gives  the  desired  results: 


f a"xcos<nx»>  = 


l an*sin(nx<p) 
n=0 


l-2a  cos(x4>)+a 

+axsin(x<t>) 

l-2axcos(x<(>)+a2x 


Also,  from  A1 , 


v anx  _ a 

^ x 

n=l  1-a 


Integrating  from  -»  to  x gives 


I pp  = - £rt(l-aX)  . 
n=l  n 

Replacing  a by  aeJ*  and  equating  real  parts  gives 


“ .nx  1 y ?y 

l -jp  cos(nxij))  = - j fM[l-2axcos(x<f>)+a^x] 
n=l 


1 


